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Abstract. The complex reflected and transmitted amplitudes from a Fabry-Perot interferom¬ 
eter are analyzed using a phase-space approach, in which the real and imaginary parts of those 
amplitudes are taken as basic variables. As functions of the optical path length of the cavity, 
the reflected amplitude describes a circle, whereas the transmitted is given by a hippopede. 
The system performance can be directly related to the geometrical parameters of these curves. 


1. Introduction 

The Fabry-Perot (FP) interferometer provides a superb illustration of the mysterious ways 
in which interference works. Despite its apparent simplicity, it plays a central role as high- 
resolution spectrometer, laser resonator, or spectral filter, to cite only but a few of its many 
relevant uses. A complete account of the subject can be found in the two comprehensive 
monographs by Hernandez |[T) and Vaughan 0. 

This variety of fields of application spawned many descriptions of the FP operation, 
each one capitalizing on specific aspects. The geometric treatment, in which one adds the 
multiple beams reflected at each of the different interfaces, is probably more instructive and, 
accordingly, is reproduced in almost every textbook Q. The question can also be tackled 
by imposing the appropriate boundary conditions, which gives the resonant frequencies and 
allowed fields in the FP Q. As the boundary conditions appear as a linear system, they can 
pop up under multiple guises; for example, the FP may be viewed as an optical transmission 
system with feedback Q, or as a direct application of the transfer matrix ||^|^, a method 
especially germane to deal with layered structures. 

Irrespective of the approach, the focus is always on the intensity distributions (both 
in reflection and transmission); namely, the well-known Airy formulas. Even if the 
corresponding amplitudes are somehow required to determine these distributions, they are 
ultimately discarded on the grounds that real experiments measure intensity. 

This viewpoint can be challenged, however, by a discussion of the conventional harmonic 
oscillator, wherein the complex amplitude is decomposed into two orthogonal in-phase and 
out-of-phase quadratures, also known as the dispersive and absorptive components Q. These 
quadratures are of paramount importance as they convey more useful information than just 
the intensity. At the quantum level, for example, they play the role of the effective position 
and momentum of the oscillator Q. Actually, a good deal of the latest advances in quantum 
information processing stem from a proper engineering of these quadratures, with homodyne 
detection constituting an ideal tool for their measurement, whereas squeezing them provides 
an efficient route to producing entanglement ID- 

Inspired by this, we intend to shed light on the amplitude response of the FP. Indeed, 
one can define an equivalent version of the quadratures. When the parameters of the FP 
vary, the amplitudes trace out elementary curves: the reflected amplitude is a circle, and the 
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transmitted one is a hippopede, a curve with remarkable properties |n][Tg. Furthermore, the 
FP performance can be naturally assigned to the geometrical features of these curves. 

Apart from its elegance, this approach makes a close contact with phase-space methods 
that pervade physics today. The derivation is straight and easy; suitable for undergraduates. 
Surprisingly enough, to the best of our knowledge, such simple ideas have been not hitherto 
explored; even if they can be important in some instances in which the phase introduced by 
the FP matters, as it happens in optical metrology, where the stabilization is crucial. 


2. The Fabry-Perot: Basic background 


The ideal FP interferometer consists of two parallel mirrors (that, for simplicity, we assume 
to be identical) separated at a distance d. Figure [T] is a block diagram of the system. This 
can be addressed by considering a plane parallel plate of thickness d and refractive index 
n immersed in a medium of index n'. The plate is illuminated near normal incidence with a 
linearly polarized quasi-monochromatic plane wave, with the electric field lying either parallel 
or perpendicular to the plane of incidence. Any diffraction effect or polarization dependence 
are thus neglected. 

The complex reflection and transmission coefficients (i.e., the ratios of the reflected and 
the transmitted amplitudes to the incident one, respectively) are given by |[^ 


Rm = 


r[l — exp(— 

1 — r^ exp{—i2^) 


T{^) = 


(1 — r^)exp(—/4>) 
1 — r^ exp (—12^) 


( 2 . 1 ) 


Here, r is the Fresnel reflection coefficient for a wave travelling from the surrounding medium 
into the FP and 


2n 

•P = -^ndcos 9 

A 


( 2 . 2 ) 


is the plate phase thickness, with A the wavelength in vacuum and 9 the angle of refraction in 
the medium n, which is related to the angle of incidence according to Snell’s law. 

The usual analysis proceeds by calculating the reflectivity and transmissivity (i.e., the 
ratios of the reflected and the transmitted intensities to the incident one). The expressions are 



Figure 1. Schematic block diagram of an FP. The intracavity fields can be analyzed from a 
variety of perspectives. The interferometer plate surfaces have reflection coefficient r and they 
are separated at a distance d. The reflected and transmitted amplitudes are labeled by R and T, 
respectively, since we assume unit incident amplitude. 
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Figure 2. Transmissivity 3^ of the FP as a function of the phase shift <I> and the parameter |r| 
whose square is the reflectivity of the plate surfaces. 


obtained directly from equation and read 

^ , ,9 Fsin^O 

fg=\R\^ = 


1 + F sin <I> 




lrp = 


1 


1 +Fsin O 


where the parameter F is 


F = 


4k|- 


(l-|r|2)2 


(2.3) 


(2.4) 


Although r is up to now a real number, we formally treat it as a complex for reasons that will 
become apparent soon and so |r|2 is the reflectivity of the plate surfaces. Equations (2.31 
constitute the time honored Airy formulas. Evidently, since there are no losses, the two 
patterns are complementary, in the sense that 


= (2.5) 

In figure 1^ we plot the transmissivity ^ as a. function of the phase thickness <I> and |r|. 
As |r| increases, the minima of fall and the maxima become sharper. In the limit of high 
|r|, the pattern consists on narrow bright fringes on an almost completely dark background. 

The sharpness of the fringes is conveniently measured by their full width at half 
maximum (EWHM), which is the width between the points on either side of a maximum 
where the intensity has fallen to half its maximum value. The ratio of the separation of 
adjacent fringes (also called free spectral range) and the EWHM is called the finesse of 
the fringes. A direct calculation shows that 




ttv/f 

2 


( 2 . 6 ) 


This quantity is a measure of the apparatus ability to resolve closely spaced spectral features. 
High values of ^ require an increased reflectivity jrp and this is accomplished by coating the 
plates surfaces with a mirror. In what follows, we assume that such a mirror is lossless. In 
that case, the Airy formulas still hold provided we interpret r as the reflection coefficient of 
the mirror (which becomes now a complex number). This adds to the plate phase thickness 
<I> a phase change on the reflection at the mirrors. In general, both modulus and phase of 
the complex r depend on the angle of incidence and the dispersion properties of the material, 
albeit such a variation can be disregarded for most practical purposes. 








Amplitude response of a Fabry-Perot interferometer 

3. Amplitude response of the Fabry-Perot 


4 


Let us look in more detail at the amplitudes (2.1 1 . First of all, we observe that is a K- 
periodic function, while T ('!>) is 27r-periodic. Such a difference cannot be seen in the intensity 
response, for both and ^ have the same period n. 

To proceed further, let us rewrite equation as 


, , 2|r|sin^r / x / x-i 

jv(o) + 


r (o) = [x{^)-iY{<t>)]. 


(3.1) 


/A(4>) 


where we have defined 


X(<I>) = (1 + |r|2)sin4>, y(4>) = (1 - |r|2)cos<I>, 

N{^) = (1 - |r|2)2+4|rpsin2 4>. 


(3.2) 


We recall that any harmonic signal .r(f) can be decomposed as 

x(f) = Acostut+ y sinojf, (3.3) 


where the X and Y are the in-phase and (7t/2) out-of-phase quadratures. In this spirit, A(<I>) 
and y(4>) can be seen as sort of quadratures for /?(•!>) and r(0). We stress, however, that 
at difference of the harmonic oscillator, here the pre-factors in their dehnition p.l| i are not 
constant, but depend on <I>. This arises because the complex amplitudes Rfi>) and do 
not have constant modulus, as in the oscillator. 

In figure[^we represent these complex amplitudes for several values of |r|. The different 
behavior commented above in relation with the periodicity translates into the fact that when 
completes a revolution, makes two turns. Both amplitudes lie inside the unit disk 
because of the fundamental constraint 




Figure 3. Phase-space trajectories for^(<J>) (left) and (right) as given in equation 
They lie inside the unit disk. The different curves correspond to different values of the plate 
reflectivity ranging from r = 0.11 to r = 0.99 in steps of 0.11 For R{^), the curves increase in 
size with |r|, while the converse happens for T(<I>). 
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To get a better grasp of these amplitudes, we introduce polar coordinates as 
|/;|exp(/p), r = |7’|exp(/T). 

With this parametrization, we can recast equations as 
,^ = 4a^cos^p, = 1 — sin^ T, 


where we have used the definition in (2.3 1 and a is a real parameter 


i + kl- 


(3.4) 

(3.5) 

(3.6) 


so that 0 < a< 1 / 2 . 

In this way, the reflected amplitude can be immediately identihed as a circle of 
radius a centered in the point {a, 0 ) of the real axis. 

The transmitted amplitude 7’(<I>) describes a hippopede (which literally means “horse 
fetter”). It was hrst investigated by Proclus HD and later on by Booth p4| , hence their 
names are sometimes attached to this stunning curve. For 0 < a < 1 /s/% it is an oval, and 
for 1 / 1/8 < a < I /2 it is an indented oval, which tends to be an eight in the limit a = I/2 
(i.e., |r| —>^ I). Yet not so well known in Physics, it has a truly amazing set of properties that 
the reader can look up in the abundant literature on the subject We merely quote 

that the hippopede can be defined as the curve formed by the intersection of a torus and a 
plane parallel to the axis of the torus and tangent to it on the interior circle. It is thus a spiric 
section m- 

For every value of |r|, the reflected amplitude passes through the origin: is zero for 

<I> = 0 and 7Z and traces the ci rcle c lockwise, getting its maximum at <I> = ±;r/2, where p = 0 
and then, according equation (3.5 i, fimax = 4a^. 

On the other hand, the transmitted amplitude also describes the hippopede clockwise. 
At <I> = 0 and n, 7’(<I>) reaches its maxima, which are in the real axis at T = I and —I, 
respectively. The minimum occurs at '!> = n/l and 37r/2, where T = —n/l and —hn/l, 
respectively. Therefore ,5jnin = I ~ which corresponds to half the waist of the hipoppede 
at its indentation. 



Figure 4. Reflected (red) and transmitted (blue) amplitudes R(0) and ^(O) for a transparent 
FP with I r\ — 0.77. For every fixed value of the parameter <I> the corresponding position vectors 
are orthogonal. 
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Figure 5. Transmitted phase lag z by an FP as a function of <I> for different values of |r|, from 
0.11 (yellow full line) to 0.99 (black full line) in steps of 0.22. The curves bend more as |r| 
increases. 


One can check that 
Ti<t>) 


= iy/p sinO, 


(3.7) 


which, in turn, implies that for a transparent symmetric system, as the one we are dealing 
with, we have 

n 


P(0)-t(0) = ±- 


(3.8) 


Consequently, for every value of <I> the reflected and transmitted amplitudes are at quadrature. 
This is illustrated in figure where we see that the position vectors of R{^) and 7’(<I>) are 
orthogonal, as it is somehow implicit in equation p.l| l. Apropos of this, it is noteworthy 
to mention that quite similar relations may be derived under the general assumptions of 
symmetric and lossless systems 119 .|2T[. 

We next examine the local slopes p(4>) and the dot denoting derivative respect to 
the parameter. They are the “rates” at which the curves are traced out. Indeed, they entail a 
valuable physical interpretation. If we focus for simplicity at T, we can write 


(3.9) 


dT 

with 0) being the angular frequency. Now, d^jdat is the single-pass time inside the cavity 
medium (for a non-dispersive material, this is nd cos 0/c, where c is the velocity of light 
in vacuum), and dz/d(0 is the time flight through the FP, which incorporates the feedback. 
Hence, dx/d(0 may be viewed as an enhancement factor of the time flight due to the 
FP ||22]|23). 

Because of equation (3.8i, both are equal for a lossless medium: p(4>) = In 

addition, we have 


T 


1 + 1 


1 - 


.12 ’ 


(3.10) 


where the negative sign indicates that the curve is oriented clockwise. This quotient is thus 
independent of <I>: where the transmitted amplitude is large, so is the velocity and the opposite. 
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Figure 6. Same as in figure]^ but for an absorbing FP made of germanium with complex 
refractive index = 5.588 — i 0.933 at a wavelength of 0.6199 ^um. We take normal incidence 
and the film thickness d varying between 0 and 0.35 /im. The marked position vectors 
correspond iod — 0.052 ^m. 


To gain further insight into this issue, in figurej^we plot T as a function of <I> for several 
values of |r|. Note that the range of variation of T is from —2n to 0 when <I> goes from 0 
to 2k, but p ranges from — ;r/2 to ;r/2, as one can directly infer at a glance from figure]^ 
When |r| is small, T is almost a straight line, with a slope pretty constant; this is the case 
when the hippopede is almost an oval, without indentation. However, as |r| increases, T starts 
bending near O = —7r/2, which is precisely at the waist. In the limit |r| —>^ 1, T becomes 
almost horizontal, with slope zero almost everywhere, except a narrow interval around the 
maxima of where it quickly gets large. 

The previous considerations can be extended to a lossy (or gain) cavity medium, a case 
of particular interest in laser physics. The medium is now specified by a complex refractive 
index (and so a complex <I>), whose imaginary part accounts for the losses. The resulting 
trajectories, for the simple case of a plate of germanium, are shown in figure]^ and they turn 
out to be loxodromics p4| , a universal feature of absorption. They start at = 0 and T = \ 
(when there is no film) and tend to R ^ r and T ^ 0 (when the film becomes opaque). We 
can appreciate that the orthogonality between position vectors does not hold true, neither the 
complementary relation equation p.5| l. 

4. Concluding remarks 

In summary, we have thoroughly explored the amplitude response of an FP interferometer. 
Despite its basic nature and its simplicity, this topic has been snubbed in the literature, which 
has strengthened instead the role of the associated intensity. 

Given the relevant role played by the Airy formulas, one might have expected that their 
“square root” counterparts should have remarkable features. This is indeed the case, as 
our results indicate: the reflected amplitude traces a circle, while the transmitted one is a 
hippopede, an intriguing curve full of nice mathematical properties. 

We finally stress that the phase-space methods employed here are quite appealing for 
they have branched into offshoots of importance for modern physical theories. 
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